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^D • Abstract. The first part of this report reviews recent developments at the interface between lattice 

^N ' work on QCD with light dynamical quarks, effective field theory and low energy precision exper- 

^H . iments. Then I discuss how dispersion theory can be used to analyze the low energy structure of 

the nn scattering amplitude in a model independent manner. This leads to an exact formula for the 



Oh. 



■^ \ mass and width of the lowest few resonances, in terms of observable quantities. As an application, I 
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consider the pole position of the d, paying particular to error propagation in the numerical analysis. 
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CSJ \ QCD with massless quarks is the ideal of a theory: it does not contain a single dimen- 

QQ ■ sionless free parameter. At high energies, the degrees of freedom occurring in the La- 

grangian are suitable for a description of the phenomena and the interaction among these 
degrees of freedom can be treated as a perturbation. At low energies, on the other hand, 

(^ ■ QCD reveals a rich spectrum of hadrons, the understanding of which is beyond the reach 

OO . of perturbation theory. In my opinion, one of the main challenges within the Standard 

Model is to understand how an intrinsically simple beauty like QCD can give rise to the 
amazing structures observed at low energy. 

^ ■ The progress achieved in understanding the low energy properties of QCD has been 

very slow. A large fraction of the papers written in this field does not concern QCD 
as such, but models that resemble it in one way or the other: constituent quarks, NJL- 
model, linear a model, hidden local symmetry, AdS/CFT and many others. Some of 
these may be viewed as simplified versions of QCD that do catch some of the salient 
features of the theory at the semi-quantitative level, but none provides a basis for a 
coherent approximation scheme that would allow us, in principle, to solve QCD. 

This talk concerns the model independent approach to the problem based on disper- 
sion theory. More precisely, I would like to show that the low energy structure of the 
sector with the quantum numbers of the vacuum, I = i = 0, can quantitatively be under- 
stood on the basis of the symmetries of QCD, without invoking any model, but instead 
relying on the available, rather crude experimental information at energies above 1 GeV. 
At low energies, the main characteristic of QCD is that the energy gap is very small, 
Myi ~ 140 MeV. More than 10 years before the discovery of QCD, Nambu [2] found 
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out why that is so: the gap is small because the strong interactions have an approximate 
chiral symmetry. Indeed, QCD does have this property: for yet unknown reasons, two of 
the quarks happen to be very light. The symmetry is not perfect, but nearly so: m^ and 
Md are tiny. The mass gap is small because the symmetry is "hidden" or "spontaneously 
broken": for dynamical reasons, the ground state of the theory is not invariant under 
chiral rotations, not even approximately. The spontaneuous breakdown of an exact Lie 
group symmetry gives rise to strictly massless particles, "Goldstone bosons". In QCD, 
the pions play this role: they would be strictly massless if m„ and m^ were zero, so 
that the symmetry would be exact. The only term in the Lagrangian of QCD that is 
not invariant under the group SU(2)xSU(2) of chiral rotations is the mass term of the 
two lightest quarks, muuu + niddd. This term equips the pions with a mass. Although 
the theoretical understanding of the ground state is still poor, we do have very strong 
indirect evidence that Nambu's conjecture is right - we know why the energy gap of 
QCD is small. 

LATTICE RESULTS RELEVANT FOR LOW ENERGY QCD 

As pointed out by Gell-Mann, Oakes and Renner [3], the square of the pion mass is 
proportional to the strength of the symmetry breaking, M^ ©« (m„ H-mj). This property 
can now be checked on the lattice, where - in principle - the quark masses can be varied 
at will. In view of the fact that in these calculations, the quarks are treated dynamically, 
the quality of the data is impressive. The masses are sufficiently light for ;fPT to allow 
a meaningful extrapolation to the quark mass values of physical interest. The results 
indicate that the ratio M^/{mu + md) is nearly constant out to values of ma^mj that are 
about an order of magnitude larger than in nature. 

The Gell-Mann-Oakes-Renner relation corresponds to the leading term in the expan- 
sion in powers of the quark masses. At next-to-leading order, this expansion contains a 
logarithm: 

l+^^^^2^1n^+0(M4)j , M^^B{m, + ma) . (1) 

Chiral symmetry fixes the coefficient of the logarithm in terms of the pion decay constant 
Fj^, but does not determine the scale A3. A crude estimate was obtained more than 20 
years ago [4], on the basis of the SU(3) mass formulae for the pseudoscalar octet. The 
result is indicated at the bottom of the left panel in Fig.l. The other entries represent 
recent lattice results for this quantity [5]-[10], which are considerably more accurate. 
The right panel shows the results for the coupling constant £4, which determines the 
quark mass dependence of the pion decay constant at NLO of the chiral expansion. In 
that case, we obtained a rather accurate result in 2001, from a dispersive analysis of the 
scalar pion form factor [11]. The lattice determinations of £4 have reached comparable 
accuracy and corroborate the outcome of our analysis. 

The hidden symmetry not only controls the size of the energy gap, but also determines 
the interaction of the Goldstone bosons at low energies, among themselves, as well as 
with other hadrons. In particular, as pointed out by Weinberg [12], the leading term in the 
chiral expansion of the S-wave nn scattering lengths (tree level of the effective theory) 
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FIGURE 1. Determinations of the effective coupling constants £3 and Ha,- 



is determined by the pion decay constant. The corresponding numerical values of «[] 
and Oq are indicated by the leftmost dot in Fig. 2, while the other two show the result 
obtained at NLO and NNLO of the chiral expansion, respectively. The exotic scattering 
length Oq is barely affected by the higher order corrections, but the shift seen in «[] is 
quite substantial. The physics behind this enhancement of the perturbations generated 
by m„ and m^/ is well understood: it is a consequence of the final state interaction, which 
is attractive in the 5*^-wave, rapidly grows with the energy and hence produces chiral 
logarithms with large coefficients. 

Near the center of the Mandelstam triangle, the contributions from higher orders 
of the chiral expansion are small [11]. Using dispersion theory to reach the physical 
region, we arrived at the remarkably sharp predictions for the two scattering lengths 
indicated on the left panel of Fig. 2. Our analysis also shows that the corrections 
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FIGURE 2. Comparing the theoretical predictions for the nn S-wave scattering lengths with lattice 
results (left) and with experiment (right). 
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to Weinberg's low energy theorem for 0.^,0^ are dominated by the effective coupling 
constants £3,14 discussed above - if these are known, the scattering lengths can be 
calculated within small uncertainties. Except for the horizontal band, which represents 
a direct determination of Uq based on the volume dependence of the levels [13], all of 
the lattice results for the scattering lengths shown on the left panel of Fig. 2 are obtained 
in this way from the corresponding results for £3 and £4. The figure neatly demonstrates 
that the lattice results confirm the predictions for Uq^Uq. 

PRECISION EXPERIMENTS AT LOW ENERGY 

The right panel of Fig. 2 compares the predictions for the scattering lengths with recent 
experimental results. While the Ke4 data of E865 [14], the DIRAC experiment[15] and 
the NA48 data on the cusp in ^ — * 3;r [16] all confirm the theoretical expectations, the 
most precise source of information, the beautiful ^^4 data of NA48 [17], gives rise to 
a puzzle. The Watson theorem implies that - if the electromagnetic interaction and the 
difference between niu and nid are neglected - the relative phase of the form factors 
describing the decay K -^ evnn coincides with the difference 5q — 5/ of scattering 
phase shifts. At the precision achieved, the data on the form factor phase do not agree 
with the theoretical prediction for the phase shifts. 

The origin of the discrepancy was identified by Colangelo, Gasser and Rusetsky [18]. 
The problem has to do with the fact that neutral kaons may first decay into a pair of 
neutral pions, which then undergoes scattering and winds up as a charged pair. The mass 
difference between the charged and neutral pions affects this process in a pronounced 
manner: it pushes the form factor phase up by about half a degree - an isospin breaking 
effect, due almost exclusively to the electromagnetic interaction. Fig. 3 shows that the 
discrepancy disappears if the NA48 data on the relative phase of the form factors are 
corrected for isospin breaking. Accordingly, the range of scattering lengths allowed by 
these data, shown on the right panel of Fig. 2, is in perfect agreement with the prediction. 
The intersection of this range with the band from the low energy theorem for the scalar 
radius (left panel) yields a'^ = 0.220(9). 

In the mass range Mj^: > 350 MeV, Fig. 3 indicates a marginal disagreement between 
NA48/2 and E865. While E865 collects all events in this region in a single bin, the 
resolution of NA48/2 is better. The fit to all Kg4 data is therefore dominated by NA48/2. 
For a detailed discussion of these issues, I refer to the talks by B. Bloch-Devaux, G. 
Colangelo and J. Gasser at KAON 2007. 1 conclude that the puzzle is gone: Ke4 confirms 
the theory to remarkable precision. 

ROY EQUATIONS 

As mentioned above, the straightforward evaluation of the chiral perturbation series 
for the nn scattering amplitude is useful only in a very limited range of the kinematic 
variables - definitely, the resonance poles are outside this range. The domain of validity 
can be extended considerably by means of dispersion theory: analyticity, unitarity and 
crossing symmetry essentially determine the low energy properties of the scattering 
amplitude in terms of the two S-wave scattering lengths. I now wish to show that the 
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FIGURE 3. Comparison of Kgn, data with the prediction for 8^ — 5/ 



properties of the lowest resonance of QCD can be worked out on this basis. 

From the point of view of dispersion theory, tztz scattering is particularly simple: 
the S-, t- and M-channels represent the same physical process. As a consequence, the 
scattering amplitude can be represented as a dispersion integral over the imaginary part 
and the integral exclusively extends over the physical region [19]. The projection of 
the amplitude on the partial waves leads to a dispersive representation for these, the 
Roy equations. I denote the S-matrix elements by S[ = rj'^explidj^ and use the standard 



normalization for the corresponding partial wave amplitudes t^: 



S',{s) = l+2ip{s)4{s), p{s) = Jl-AMl/s. 



(2) 



The S-matrix elements and the partial wave amplitudes are analytic in the cut 5-plane. 
There is a right hand cut (4M^ < i' < oo) as well as a left hand cut (—0° < s < 0). The 
Roy equation for the partial wave amplitdue with 1 = 1 = reads 

t^^s)=a + {s-ml)b+Y.I. ds'K'o'i^,s')lm4is'). (3) 

The equation contains two subtraction constants. As a consequence of crossing symme- 
try, these can be expressed in terms of the S-wave scattering lengths: 



a = a, 



0' 



Z7=(2a[j-5a^)/12M2. 



(4) 



The kernels K^/^, {s, s') are explicitly known algebraic expressions which only involve the 
variables s,s' and the mass of the pion. The integrals on the right hand side thus only 
involve observable quantities: the imaginary parts of the partial waves. 

As demonstrated by Roy, his equations rigorously follow from general principles of 
quantum field theory. They are valid for real values of 5 in the interval — 4M^ < 5 < 60A/^ 
(the upper end is pushed up to 68M^ if the scattering amplitude obeys Mandelstam 
analyticity). Using known results of general quantum field theory [20, 21], we have 
shown that these equations also hold for complex values of s, in the intersection of the 
relevant Lehmann-Martin ellipses [1]. 
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The pioneering work on the physics of the Roy equations was carried out more than 30 
years ago [22]. The main problem encountered at that time was that the two subtraction 
constants were not known. These dominate the dispersive representation at low energies, 
but since the data available at the time were consistent with a very broad range of 5- wave 
scattering lengths, the Roy equation analysis was not conclusive. 

The insights gained by means of ;fPT thoroughly changed the situation. As discussed 
in detail in the first part of this report, the S-wave scattering lengths are now known very 
accurately. The main limitation in the numerical evaluation of equation (3) arises from 
the accuracy to which the imaginary parts can be pinned down. In this connection, it is 
essential that the Roy equations involve two subtractions, so that the kernels fall off with 
the third power of the variable of integration. This ensures that the contributions from 
the low energy region dominate. Note that the left hand cut plays an important role here: 
taken by itself, the part of the kernel that accounts for the right hand cut falls off only 
with the first power of the variable of integration, but the high energy tail is cancelled by 
the contribution from the left hand cut. 

At low energies, the S- and P-waves dominate. In [23], we solved the Roy equations 
for these waves only below 800 MeV, relying on the literature to estimate the contri- 
butions from higher energies and from the partial waves with i >2. In the meantime, 
we have extended our analysis and now solve the Roy equations on their full range of 
validity, not only for the S- and P-waves, but also for the D- and F- waves. Treating the 
S-wave scattering lengths, the imaginary parts above Smax = 68M^ and the elasticities 
as input, the Roy equations admit a two-parameter family of solutions. We identify the 
two free parameters with the values of the phase shifts 5q and 5/ at 800 MeV. In view 
of the excellent experimental information about the vector form factor of the pion, the 
value of 5/ (800 MeV) is reliably known, but the phenomenological information about 
5q (800 MeV) is comparatively meagre - this currently represents the main source of 
uncertainty in low energy nn scattering and will be discussed in detail below. 

POLE FORMULA 

The positions of the poles represent universal properties of QCD, which are unambigu- 
ous even if the width of the resonance turns out to be large, but they concern the non- 
perturbative domain, where an analysis in terms of the local degrees of freedom - quarks 
and gluons - is not in sight. One of the reasons why the values for the pole position of 
the o quoted by the Particle Data Group cover a very broad range is that all but one of 
these either rely on models or on the extrapolation of simple parametrizations: the data 
are represented in terms of suitable functions on the real axis and the position of the 
pole is determined by continuing this representation into the complex plane. If the width 
of the resonance is small, the ambiguities inherent in the choice of the parametrization 
do not significantly affect the result, but the width of the o is not small. For a thorough 
discussion of the sensitivity of the pole position to the freedom inherent in the choice of 
the parametrization, I refer to [24]. 

The determination of the o pole provides a good illustration for the strength of the 
dispersive method and for the relative importance of the various terms on the right hand 
side of the Roy equations. The representation of the S-matrix element given above holds 
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on a limited region of the first sheet. The pole sits on the second sheet, which is reached 
from the first by analytic continuation from the upper half plane into the lower half plane, 
crossing the real axis in the interval 4M^ < s < 16M^, where the scattering is elastic. 
Now, unitarity implies that the values of the S-matrix element on the first and second 
sheets are related by 5'q(5)^^ = l/S^{sy . Hence a pole on the second sheet occurs if and 
only if 5q(5) has a zero on the first sheet. Accordingly, we have an exact equation, which 
allows us to study the behaviour of the amplitude in the vicinity of the pole and find out 
whether - in the limited region of the first sheet where the Roy equations are valid - 
there are any resonances with the quantum numbers of the vacuum: 

S^o{s) = 0. (5) 

As emphasized above, the representation of the S-matrix element that follows from 
the Roy equations exclusively involves observable quantities and can be evaluated for 
complex values of s just as well as for real values - for the above formula, an analytic 
continuation is not needed. 

Inserting our central representation for the scattering amplitude in (5), we find that, 
in the region where the Roy equations are valid, the function Sq{s) has two zeros in the 
lower half of the first sheet: one at ^/s = 441 — z 272 MeV, the other in the vicinity of 1 
GeV [1]. While the first corresponds to the state /o(600), commonly referred to as the o, 
the second zero represents the well-established resonance /o(980). Our analysis sheds 
little light on the properties of the latter, because the location of the zero is sensitive to 
the input used for the elasticity r]^ (5) - the shape of the dip in r^Q (5) and the position of 
the zero represent two sides of the same coin. For this reason, I only discuss the (7. 

DISCUSSION 

We are by no means the first to find a resonance in the vicinity of the above position. In 
the list of papers quoted by the Particle Data Group [25], the earliest one with a pole in 
this ball park appeared more than 20 years ago [26]. What is new is that we can perform 
a controlled error calculation, because our method is free of the systematic theoretical 
errors inherent in models and parametrizations. For this purpose, it is convenient to split 
the right hand side of the Roy equation for ?q (5) into three parts: 

1. Subtraction terms 

2. Contribution from Im?Q (5) below KK threshold 

3. Contributions from higher energies and other partial waves 

ad 1. As discussed above, the subtraction terms are determined by the S-wave scatter- 
ing lengths. The prediction of [11] reads a^ = 0.220 ±0.005, a^ = -0.0444 ±0.0010. 
Following the propagation of a change in a^, we find that an increase by 0.005 shifts the 
pole position by (—2.4 ± z 3.8) MeV, while the response to an increase in Uq by 0.0010 is 
a shift of (0.8 — /4.0) MeV [1]. These numbers show that the error in the pole position 
due to the uncertainties in the subtraction constants are small. 

ad 2. Below KK threshold, the S-waves are elastic to a very good approximation. The 
function Im?Q (5) shows a broad bump, nearly hits the unitarity limit somewhere between 
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FIGURE 4. Behaviour of &? below KK threshold 



800 and 900 MeV and then rapidly drops, because the phase steeply rises, reaching 180° 
in the vicinity of 2 M^. Hence there is a pronounced dip in lmtQ{s) near KK threshold: 
the behaviour of the imaginary part is controlled almost entirely by the phase shift 8q{s) . 
Fig. 4 shows several recent representations for this phase. Replacing the integral over our 
central representation for Im?Q (i') from 4M^ to 4M| by the one of Bugg [27] and leaving 
everything else as it is, the pole moves to 444 - i 267 MeV. Repeating the exercise 
with the representations of Achasov and Kiseiev [28], Albaiadejo and Oiler [29] and 
Kamihski, Pelaez and Yndurain [30], the pole is shifted to 438 - i 274 MeV, 451 - i 257 
MeV and 458 - i 253 MeV, respectively. As mentioned above, the solution of the Roy 
equation for tQ{s) depends on the value of 5a = 5q (800 MeV). The shaded band in the 
figure shows the behaviour of the Roy solutions for 78.3° < 5a < 92.3°, the range used 
in [1]. If the imaginary part of t^is) is evaluated with the lower edge of this band, the 
pole occurs at 435 - i 276 MeV, while the upper edge corresponds to 456 - i 262 MeV. 

The representation used by Yndurain and collaborators underwent a sequence of 
gradual improvements. The most recent edition is definitely better than the earlier 
versions, but an important difference to our representation remains: as can be seen on the 
right panel of Fig. 4, the phase 5q{s) of KPY III [30] still contains a kink (discontinuity 
in the first derivative) at 932 MeV, as well as a hump below that energy. The kink 
arises because (i) two different parametrizations are used below and above 932 MeV 
and (ii) the correlations between the two regions imposed by causality are ignored. The 
deficiency is discussed in detail in [3 1], where it is shown that both the hump and the kink 
are artefacts, produced by the use of a parametrization that is not flexible enough. The 
analysis in [24] fully confirms this conclusion: among the 42 conformal parametrizations 
constructed there, all of those with an acceptable behaviour above 900 MeV run within 
the band of Roy solutions specified above. 

The hump is also responsible for the disagreement between [11] and [30] regarding 
the scattering lengths of the partial waves with angular momentum £= 1,2,3 and the ef- 
fective ranges with £ = 0, 1 , 2, 3: the contributions from Im?Q (5) to the sum rules for these 
quantities are not the same (the sum rules are listed in Eqs. (14.1) and (14.3) of [23]). 
Within errors, the difference between these contributions reproduces the difference in 
the quoted results, without exception. 
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I add a remark concerning the model of Hanhart, Pelaez and Rios [32], who apply 
the inverse amplitude method to improve the one loop approximation to the chiral 
perturbation series of SU(2)xSU(2). In the original formulation of the model, the chiral 
expansion ?Q (5) = t2{s) + t4{s) + . . . is unitarizedwith?Q(5) =t2{s)/{\ —t4{s)/t2{s)},hut 
this recipe fails in the vicinity of the Adler zero, because the term ^4(5) does not vanish 
there. The deficiency is readily cured. It suffices to replace the lAM formula with 

^o(*) = -rVrT777T' h{s)=t2{s)-t2{sA,), k{s)=t4{s)+t2{sA,), (6) 

\-t4{s)/t2{s) 

where sa^ is the position of the Adler zero in one loop approximation. Since ^2(^4) 
represents a term of 0{p^), the chiral expansion of (6) reproduces the one loop approx- 
imation of ;tPT, also in the vicinity of the Adler zero. A similar recipe is used in [32]. 

The model exclusively involves the coupling constants F^.ii, ... ,£4 of the effective 
Lagrangian. As discussed above, £3 and £4 are known quite well; i[ and £2 can be 
determined on phenomenological grounds [11]. The result for the phase shift obtained 
by inserting the numerical values in the above formula is indicated on the right panel of 
Fig. 4. This shows that the model yields a decent approximation only below 500 MeV. 
The parametrization used by Hanhart eta al. [32] is better, because these authors treat the 
coupling constants ii and £2 as free parameters. This extends the range of energies where 
the lAM parametrization makes sense, but since the model does not account for the sharp 
increase in the phase towards KK threshold, it can at best give a semi-quantitative picture 
of the o. For the parameter values adopted in [32], the zero of the denominator in (6) 
occurs at 444(6) - i 218(10) MeV: the mass is OK, but the width is too low by 100 MeV. 
Inserting the observed values of ii and £2, the zero moves to 413(12) - i 269(12) MeV: 
now the width is OK, but the mass is too low. 

ad 3. Finally, I turn to the contributions of the third category: higher energies and 
other partial waves. Among these, the one from the P-wave, for example, is by no means 
negligible, but, as mentioned above, this wave is known very well. In fact, in the vicinity 
of the zero of Sq{s), the sum of the contributions of this category can be worked out 
quite accurately. In [1], we estimated the net uncertainty in the pole position from this 
source at ± 4 ± i 6 MeV. As a check, we can simply replace our central representation 
for the contributions of category 3 by the one in [30], retaining our own representation 
only for the remainder. The operation shifts the pole position by - 0.6 - i 1.2 MeV, well 
within the estimated range. 

CONCLUSION 

Adding the errors up in square, the result for the pole position becomes [1] 

^ = 44 1 + 8^ - / 272 + 1^2.5 MeV . (7) 

The error bars account for all sources of uncertainty and are an order of magnitude 
smaller than for the crude estimate y^ = (400 - 1200) - i (250 - 500) MeV quoted by 
the Particle Data Group [25]. The dispersive representation of the S-matrix element also 
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allows us to calculate the residue of the pole occurring on the second sheet, 

tS{sf=^^ + ... (8) 

S-S(y 

Our preliminary result for the magnitude of the residue is Ir^l =0.218_ooio GeV^. 
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